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We study normal modes propagating on top of the stable uniform background in arrays of dipolar 
Bose-Einstein condensate (BEC) droplets trapped in a deep optical lattice. Both the on-site mean- 
field dynamics of the droplets and their displacement due to the repulsive dipole-dipole interactions 
(DDIs) are taken into account. Dispersion relations for two modes, viz., high- and low- frequency 
' counterparts of optical and acoustic phonon modes in condensed matter, are derived analytically 

and verified by direct simulations, for both cases of the repulsive and attractive contact interactions. 
The (counterpart of the) optical-phonon branch does not exist without the DDIs. These results are 
£Nj ' relevant in the connection to emerging experimental techniques enabling real-time imaging of the 

^ , condensate dynamics and direct experimental measurement of phonon dispersion relations in BECs. 

<D 



- I - < 



PACS numbers: 03.75.Lm; 05.45.Yv; 63.20.-e; 63.20.Dj 



I. INTRODUCTION 

The propagation of collective excitations in periodically structured media has been recognized long ago as a funda- 
mentally important topic of condensed-matter physics [l[ . More recently, it has been understood that many condensed 
media, where the intrinsic dynamics is very complex, may be effectively "simulated" by rarefied quantum gases, for 
which it is much easier to predict and observe fundamental dynamical effects [2Lj3j]. In the latter context, the role 
of the periodic structure is commonly played by optical-lattice (OL) potentials In particular, the use of atomic 
Bose-Einstein condensates (BECs) trapped in OLs makes it possible to simulate the propagation of acoustic waves in 
I ■ crystals theoretically [5[-[7| and in a direct experiment [8j. 
\ A deep OL potential splits the condensate into an array of droplets trapped in local potential wells, which are 
weakly coupled with nearest neighbors by tunneling of atoms across potential barriers. In this case, it is natural to 
approximate the global BEC wave function by a superposition of wavelets of locally trapped atoms, which resemble 
well-known Wannier modes that can be used an alternative basis instead of Bloch waves @. On the basis of this 
expansion, the reduction of the underlying full Gross-Pitaevskii equation (GPE) to its discrete version (DGPE) has 
been rigorously derived in several works [loj | and reviewed in Ref. [ll[. The derivation of the DGPE was further 
extended to the model with a relatively strong contact nonlinearity, when the onsite nonlinearity in the resulting 
discrete equation is nonpolynomial fl2| . 

A significant extension of the variety of dynamical effects in atomic BECs is provided by long-range dipole-dipole 
£> | interactions (DDIs). Basic and more sophisticated properties of dipole condensates were predicted theoretically [la ], 
which was followed by the creation of such condensates in vapors of 52 Cr [TtJ • The flexibility of this setting is buttressed 
by the possibility to vary the relative strength of the DDI and contact interactions by means of the Feshbach resonance 
affecting the scattering length for colliding chromium atoms [l8| . Results obtained in this area have been reviewed in 
Ref. [H|. 

Further progress in the experimental work has recently resulted in the creation of dipolar condensates with the 
strongest possible atomic DDIs in 164 Dy [2fJ. Also promising is the work with erbium [2l|, and with gases composed 
of molecules carrying electric dipole moments [22j ■ 

OLs provide a powerful set of tools for the studies of dipolar BECs. In particular, the stabilization and destabiliza- 
tion of the condensate, trapped in deep OLs, by DDI was demonstrated theoretically and experimentally (23[. It was 
straightforward to generalize the DGPE for this case, which leads to the discrete equation combining the local onsite 
nonlinearity with the nonlocal nonlinear interaction between the sites, in the one- O [3] and two-dimensional (2D) 
[HI settings alike. Further, we have recently demonstrated that, in case the uniform (alias continuous- wave, CW) state 
of the dipolar condensate, trapped in a deep periodic potential (i.e., the discretized state with equal amplitudes at all 
sites of the lattice), is modulationally unstable, two- and three-period density patterns emerge as energy minimizers, 
the stability area being very large for the three-period pattern (24j. Modulated-density solutions are also produced 
by the DGPE for condensates fragmented by the deep OL potential in the absence of long-range interactions, but in 
that case they are completely unstable [25(. 2D supersolid structures supported by the OL in the dipolar BEC were 
predicted too [2(| , and anisotropic DDIs may give rise to striped 2D patterns even in the absence of the lattice [27j • 



2 



Beyond the aforementioned detailed studies of static patterns, a physical problem of obvious interest is the study 
of collective excitations propagating on top of stable states supported by the long-range interactions. In this context, 
a well-known result is the prediction of the roton branch of excitations in the dipolar condensates [H[ . The subject of 
the present work is the analysis of excitation modes propagating through the dipolar BEC fragmented by the deep OL 
and, accordingly, modeled by the DGPE including the long-range DDIs. The rapid progress in experimental studies 
of condensates (in particular, the development of techniques for in-situ imaging of the condensate dynamics [29l . l30l |) 
suggests that the direct observation of such modes will become possible, hence the analysis is a subject of direct 
relevance. A related promising experimental technique, based on the method of the Bragg spectroscopy [31|, was very 
recently elaborated for directly measuring phonon dispersion relations in trapped BECs 30]. In this work, we focus 
on the most fundamental case of the normal modes excited on top of the stable CW state. Still more challenging 
situations with the underlying periodically modulated density profiles will be considered elsewhere. 

The most essential result of the analysis is that the long-range DDI between condensate droplets gives rise to a high- 
frequency branch of collective excitations, which emulates the optical-phonon mode well known in condensed-matter 
physics. This excitation propagates as a displacement wave in the array of trapped BEC droplets, which resembles 
the dynamics of the Frenkel-Kontorova model 32] . Another fundamental branch of the propagating waves represents 
low-frequency density oscillations of the droplets, which is a counterpart of acoustic phonons in condensed matter. It 
is a generalization of the branch which was recently found in the limit of immobile droplets [24j . 

Excitation modes similar to those reported here may be expected in other physical settings which feature long-range 
interactions between elements of a fragmented wave field. In particular, laser illumination of atoms may induce an 
effective gravity-like long-range attraction in BEC, which, in the combination with the OL, may give rise to a supersolid 
[33| . An example in optics is offered b y th e light field trapped in an arrayed waveguide embedded into a medium 
with the nonlocal thermal nonlinearity [34j . which gives rise to a number of patterns in one and two dimensions, as 
demonstrated theoretically |35{ and experimentally 13611 . Other possible realizations of this generic setting in optics 
may use the nonlocal nonlinearity of liquid cry stals [37j, and an effective nonlocality in semiconductor waveguides 
induced by the transport of charge carriers |38| . 



II. THE MODEL 



The on-site dynamics of the dipolar-BEC droplets trapped in local potential wells of the deep OL potential is 
governed by the DGPE, which includes the onsite cubic nonlinearity and long-range DDI. Following the notation 
adopted in Ref. [2J], the discrete equation is written as 

i^rr = -c„(/„+i - /„) - c„_i(/„_i - /„) + o |/ n | 2 /„ - r - — — -/„, (1) 

where /„ is the discrete wave function, C n are inter-site coupling coefficients proportional to the overlapping integrals 
between the Wannier modes exponentially localized at adjacent sites [10(-[13|. If the coordinate of the n-th droplet is 
taken 

+ £,n , where accounts for its displacements due to the DDIs with other droplets (see Fig. [1]) , then 
the corresponding dependence of coefficients C n in Eq. ([T]) on is 

C n = CK») exp f l+^LZ^A = ccxp (-^f^) , (2) 

c = exp (—1/1), and is a coefficient of tunneling of atoms between adjacent sites separated by distance I. 

Proceeding to nonlinear terms in Eq. (fTJ), a = — 1 and +1 correspond to the attractive or repulsive onsite contact 
interaction (abbreviated below as AC or RC, respectively), and F is the relative strength of the DDI, with respect to 
the onsite nonlinearity [13| . In the present notation, T > and T < correspond to the attractive and repulsive DDIs, 
respectively. According to its definition, T may be adjusted by means of the Feshbach resonance which affects the 
strength of the local nonlinearity, and, additionally, by fixing the orientation of the magnetic moments with respect 
to the system's axis. 

The energies of the DDI, and of the trapping of the condensate in the OL potential are 

C/dd = -r ^"^"L Uol = -U cos ( 2 ™n) , (3) 

m=£n \ Xn Xm \ n 

where Uq is the depth of the periodic potential, hence the motion of each droplet is governed by the Newton's equation, 
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where p is the density of the condensate, so that the mass and momentum of the n-th droplet are p\f n \ and 
p\fn\ d^ n /dt, as shown in Fig. Q] 
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FIG. 1: (Color online) In the top panel, circles represent the array of droplets trapped in the OL potential. The bottom panel 
displays a generic distribution of droplet masses, which are proportional to areas of the circles. 



Thus, the dynamical system following from Eqs. (JXJ) , (U|) and © takes the following form: 
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' dt 
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This system was simulated by means of the Runge-Kutta method of the sixth order. 
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III. NORMAL EXCITATION MODES 



To study collective-excitation modes in the present context, we set /„ = ((/>o+g n )e~' lflt , where <j) = y/jbj [a — 2r£(3)] 
is the amplitude of the CW state, as found in Ref. [24[ [£(3) ~ 1.2021 is the value of the Riemann's zeta- function] , 
and linearize Eqs. and ([5]) with respect to small perturbations, g n and £„: 

=[-p + 2c + 2(a- 2rC(3)) g n - c (g n+1 + g n . x ) 

q>0 q>0 



d 2 (jn 
' dt 2 



-(x) ^-MT^, 



g>0 



-24T$ £ g-6 (^_ a _ g* +q + g* _ q 
q>0 



(7) 



where q is an integer index of the summation over all sites of the lattice, the asterisk stands for the complex conjugate, 
and Ii q = g n -q + g-n+q + 9n-q + 9n+q- Difference-differential equations ([7]) are Fourier-transformed into the system of 
algebraic equations, 



A + lj -B -iF 

-B A-uj -iF 

lipfa^F ((pM^F D-lo 2 



S fe (cj) 



(8) 
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where Gfc(w) and ^^(0;) are the Fourier transforms of g n (t) and £, n (t), and 



A = -4c sin 



2 / 
su\{k,q) 



— B, B 



cos(kq) 



q>0 



9>0 q 



D 



1 /2?r 



C/o + 24r^ 



C(5)-E^ 5 cos(fcg) 

9>0 



(9) 



with C(5) w 1.0369. 

Dispersion relations for propagating waves, w = are derived from Eq. ((5J) by equating the determinant of the 

system's matrix to zero, which yields a biquadratic equation, 

- (A 2 — B 2 + D)lo 2 + D(A 2 - B 2 ) - ——{A + B) = 0. (10) 

In the parameter region of interest, which corresponds to the stable CW background, Eq. (fTO)) yields four solutions 
for the frequency, uj(k) — {oJi(k),u)2(k) — — 0Ji(k), ^(k), cj 4 (fc) = —cu 3 (k)}. In the following we display the positive- 
frequency branches, the negative ones being their specular images. Further straightforward analysis of the solutions 
produced by Eq. ([5]) demonstrates that the low- and high-frequency branches represent waves of onsite density 
oscillations and droplet displacements, respectively. The latter one may also be called a dipolar-phonon mode. 

The BEC system which does not admit the droplet displacements (in other words, the corresponding trapping 
frequency is assumed to be infinite) is modeled by the first equation in ^ with = (zjj . The corresponding 
dispersion relation amounts to the single positive real root for w 2 in the parameter region where the CW background 
is stable. It represents the density-oscillation mode in the case of frozen displacements. 

The modes in the present system may be compared to Langmuir and ion-acoustic waves in plasmas [39l |. In this 
context, the displacement mode, although it originates as the dipolar-phonon wave, is a counterpart of the high- 
frequency Langmuir wave of electron oscillations in the plasma, while the density-oscillation mode is an analog of 
the ion sound. The reason for this "switch of the roles" is that the displacement wave is actually a stream of optical 
phonons, in terms of solid state [l|, while the ion sound in the plasma is carried by acoustic phonons. Further, in 
the long- wave limit the dispersion relation for the low-frequency density mode takes precisely the form of the ui(k) 
relation for acoustic phonons, see Eq. (fT3)) below. 



IV. RESULTS AND DISCUSSION 



A. Analytical considerations 



The analysis of stability conditions for the CW background in the present system demonstrates that they are 
identical to those in the previously considered model [24| with the frozen displacement degree of freedom. For the RC 
interactions \a = +1 in Eq. (JTJ) ] , the CW-stability region extends to arbitrarily weak repulsive DDIs, T < (see Fig. 
1 in Ref. 24]), while in the case of the AC interactions (a = —1) the CW has a stability region if the repulsive DDI 
is strong enough, 

|r|>|r thr |, (n) 

with the threshold value that can be found from dispersion relation (| L0[) . For /i = 1 (the chemical potential which is 
adopted in figures displayed here), r t hr ~ —1.27. To relate these results to the experimental situation, it is relevant to 
mention that as shown in Ref. [13], a typical value |T| = 0.5, relevant to our analysis (see Fig. [2]below), corresponds 
to parameter values which are realistic in experiments with 52 Cr: the scattering length ~ 2 nm and the transverse- 
confinement radius ~ 5 mm. Actually, positive and negative values of T may be varied within broad limits by means 
of the Fcshbach resonance, as shown experimentally [171 ]. 

As mentioned above, the analytical solution of dispersion equation (fTDj) demonstrates that the stable CW background 
supports two modes of the propagating excitations, viz., the low- frequency acoustic- phonon and high-frequency optical- 
phonon modes. The two dispersion curves are displayed in Fig. [5J for both signs of the contact interaction, RC and 
AC. A noteworthy peculiarity of the latter (AC) case is that the dispersion relation has a minimum at the center of 
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FIG. 2: (Color online) The top row: dispersion curves for the low-frequency acoustic phonon mode, with repulsive (a) and 
attractive (b) signs of the contact interaction (RC and AC, respectively). The bottom row: dispersion curves for the high- 
frequency optical-phonon mode, with the RC (c) and AC (d) signs of the local nonlinearity. The curves are drawn for fi = 1, 
in the first Brillouin zone, < k < 2n. The corresponding relative values of the DDI strength, F, are indicated in the plots. 
Vertical lines in panels (a) and (c) indicate frequencies predicted by the two dispersion curves at three wavenumbers, k = 1.57, 
2.5, and 4, for which numerically computed temporal spectra of lattice perturbations, excited as per Eq. (|15[) . are displayed 
below in Fig. |4] 



the Brillouin zone, provided that the repulsive DDI is relatively weak. We stress that only the situations with the 
stable CW background are shown, therefore values of T which do not satisfy stability condition ([TT]) are not included 
in panel [2{b). 

The decrease of the largest frequency of the acoustic-phonon mode with the increase of which is observed in Fig. 
[5Ja), implies that the system becomes "less stiff', for the propagation of the density waves, with the strengthening 
of the repulsive DDI. A reason for this counter-intuitive conclusion is that the model does not directly take into 
account the DDI inside each droplet. However, this may be included via a renormalization of the onsite contact 
interaction. On the other hand, Fig. &b) demonstrates that the stiffness increases with |T| in the case of the AC 
interactions (c = —1). In either cases of the RC or AC interactions, with the increase of |T| the dispersion curves for 
the low-frequency (LF) acoustic-phonon mode asymptotically approach a limit form which can be derived from Eqs. 
© and (TO]): 



<^LF(r —> — oo) = ±4csin 



\ 



/' 



2cC(3) 



q~ 3 cos(kq). 



(12) 



<J>0 



In Figs. [2ja) and (b), the dispersion curves approach the asymptotic one from above and below in the cases of the 
RC and AC interactions, respectively. 

Naturally, these results imply that the acoustic-phonon mode can propagate in the absence of the DDI (r = 0) in 
the condensate with the RC interactions only. The addition of the DDIs modifies the corresponding dispersion curve 
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FIG. 3: (Color online) Contour plots of the density, |/ n | (a), and displacement, (b), for simultaneously excited acoustic 
and optical phonon waves, in the case of the repulsive contact and dipole-dipole interactions, a — +1 and T = —2. The 
wavenumber and amplitude of the initial excitation (|15p are k = tv/2 — 1.57 and e = 0.001. The respective frequencies of 
the propagating low- and high-frequency modes are o>lf = 1.8 (a), and ojhf = 14.9 (b), which agree with those obtained from 
dispersion relations (|10p and the temporal spectra (see Fig. 2}. 



without altering it qualitatively. Note also that the acoustic-phonon mode curve takes the form of ui 
k, with the effective sound velocity that can be found from Eqs. (fTU|) and 



vk for small 



v = 0oV2c(a-2rC(3)). 



(13) 



On the other hand, the optical-phonon mode cannot propagate in the absence of the repulsive DDI, i.e., at T = 0, 
when it degenerates into uncoupled oscillations of the droplets around local minima of the lattice potential, with 
frequency u osc — (2tt/0o) yUo/p [see the horizontal black dotted line in Fig. [5] (c)]. Further, the expansion of 
the ui(k) relation for the high-frequency (HF) optical-phonon mode in the long- wave limit (for small k 2 ) exhibits a 
quadratic dispersion, similar to that for the Langmuir waves in the plasma: 



WHF(fc) 




3$c(3) |r| 



(14) 



B. Direct simulations 



The propagation of the waves through the stable CW background was tested by direct simulations of Eqs. ([7]). The 
system was initially excited by a perturbation in the form of 

f„(t = 0) = 0o + eexp(ikn), (15) 

with small amplitude e and arbitrary wavenumber k. The simulations demonstrate a superposition of the acoustic- 
phonon and optical-phonon propagating modes, as shown in Fig. |3l 

Characteristic frequencies of both types of the propagating modes, as obtained from the direct simulations, almost 
precisely fit those predicted by the dispersion relations. This is shown in Fig. [4] by means of sample temporal spectra 
representing the propagating perturbations excited as per Eq. (fT5|) with e — 0.001 and three different wavenumbers, 
k = 1.57, 2.5, and 4, in the case of the RC interactions. For each k, the creation of the acoustic-phonon and 
optical-phonon modes is observed. The central frequencies, indicated in the plots, are very close to their counterparts 
predicted by the dispersion curves, which are marked by vertical lines in Figs. [2ja,c). 

Additional simulations (not displayed here) demonstrate that the increase of the amplitude of the initial perturbation 
(e) naturally makes the spectra wider and more complex. 
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FIG. 4: Typical examples of temporal spectra of excitations in the system with a = +1, V — — 2, fi = 1, computed at a 
particular site of the lattice: (a), (c), (e) - the acoustic-phonon modes; (b), (d), (f) - the optical-phonon modes. The peak 
values of uj accurately fit their analytical counterparts obtained from Eq. (|10p . which are marked in Fig. [2] 



V. CONCLUSIONS 



This work reports a new mode of high-frequency collective excitations propagating on top of the stable background 
in the dipolar BEC trapped in the form of the array of droplets in a deep OL potential. This branch of the excitations, 
which is similar to the optical-phonon mode in condensed matter, and also features similarity to Langmuir waves in 
plasmas, does not exist in the absence of the repulsive dipole-dipole interactions. The dispersion relations for the new 
mode and the additional low-frequency one of the acoustic-phonon type, which accounts for density oscillations in the 
trapped BEC, were derived analytically. The results have been confirmed by direct simulations of the excited system. 

The currently developin g te chniques for in situ imaging of dynamical phenomena and recording phonon dispersion 
relations in trapped BEC [29l |30| open the way to experimental realization of the results reported in this work in 
atomic gases of chromium, dysprosium and erbium, and in molecular gases, such as KRb and NaLi. As concerns the 
theoretical analysis, it can be extended in several directions. First, it would be interesting to study the nonlinear 
dynamics of excitations with nonsmall amplitudes. Second, a natural extension is to consider collective excitations 
on top of stable periodically modulated patterns, such as the double- and triple-periodic ones [24 |. This extension is 
expected to produce a larger number of different dynamical branches. Then, a challenging problem is to find normal 
modes of collective excitations in 2D settings. Finally, quite an interesting possibility is to introduce a deformable 
Bose-Hubbard model, as a quantum version of the mean-field system considered in the present work, which may be 
then implemented in chains of trapped ultracold atoms or ions [3j . 

A.M., G.G., and Lj.H. acknowledge support from the Ministry of Education and Science of Serbia (Project 0145010). 
The work of B.A.M. and T.P. was supported in a part by the German-Israel Foundation through grant No. 1-1024- 
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